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Abstract 

We study multiplicity of the supercritical traveling front solutions for 
scalar reaction-diffusion equations in infinite cylinders which invade a lin- 
early unstable equilibrium. These equations are known to possess travel- 
ing wave solutions connecting an unstable equilibrium to the closest stable 
equilibrium for all speeds exceeding a critical value. We show that these 
are, in fact, the only traveling front solutions in the considered problems 
for sufficiently large speeds. In addition, we show that other traveling 
fronts connecting to the unstable equilibrium may exist in a certain range 
of the wave speed. These results are obtained with the help of a variational 
characterization of such solutions. 

1 Introduction 

Front propagation is a ubiquitous feature of many reaction-difFusion systems 
in physics, chemistry and biology [l]-[6j. A question of particular interest that 
arises in a number of applications is related to fronts that invade an unstable 
equilibrium (for a review, see |7|). One of the well-known examples of such a 
system is given by a scalar reaction-diffusion equation in a cylindrical domain: 

ut = Au + fiu,y). (1.1) 

Here u = u{x,t) e M is the dependent variable, a; G S C M" is a point in a 
cylindrical domain S = x M, whose cross-section ft C M"^^ is a bounded 
domain with the boundary of class (not necessarily simply-connected), and 
t S [0, -l-oo) is time. We will often write x = (y, z) £ S, where y G fl is 
the coordinate on the cylinder cross-section and z G M is the coordinate along 
the cylinder axis. The function / iMxil— j'Eisa nonlinear reaction term. 
On different connected portions of the boundary we assume either Dirichlet 
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or Neumann boundary conditions (see [s] for the motivation of this particular 
choice of boundary conditions) 



I as J 



V ■ Vm 



0, 



(1.2) 



where 9E± = dVL± x ]R and 9S]o ~ x K. Furthermore, we assume that 



is an unstable solution of (1.1) in the sense of linear stability. 



By front solutions of (1.1 1, one understands traveling wave solutions, which 
are special solutions of (1.1 1 of the form u{y,z,t) = u{y,z — ct), connecting 
monotonically two distinct equilibria, i.e., stationary z-independent solutions of 
In the present context, of special interest are the fronts which invade the 
u = equilibrium. Assuming the invasion happens from left to right, the profile 
u > of such fronts is described by the following equation: 



Au + cu^ + f{u,y) ^ 0, 



lim u{- 



z) = uniformly in fl, 



(1.3) 



with the same boundary conditions as in ( |1.2[ ), where c > is the propagation 
speed. These solutions are known to play an important role in the long-time 
behavior of solutions of ( |1.1[ ) (for an overview, see e.g., [o]). 



The problem of existence and uniqueness of solutions to (1.3 1 has a long 
history, going back to the pioneering works of Kolmogorov, Petrovsky and 
Piskunov [lO|, and within the present context has been extensively studied by 
Beretsycki and Nirenberg 11 , and Vega 12 13 . Specifically, Berestycki and 



Nirenberg established existence, monotonicity, and uniqueness of solutions of 
(1.31 connecting zero to the smallest positive equilibrium for all c > cl, for 
some cl > (for the precise definition, see Sec. [2]), in the case of Neumann 
boundary conditions 



11 



Theorem 1.9]. Their methods were later extended by 
Vega to problems with Dirichlet boundary conditions 12 Theorem 3.4(B)]. We 



note that in these works uniqueness of the front solutions is understood in the 
class of functions connecting zero ahead of the front to a prescribed equilibrium 



behind the front. In general, front solutions of (1.3) may fail to be unique, even 
for a fixed value of c, if they are allowed to approach different equilibria. One 



example of multiplicity of the front solutions to (1.3 1 was given by Vega under 



specific assumptions on the set of equilibria of (1.1) 12, Remark 3.5(B)]. 

The question of multiplicity of solutions of (1.3 1 and (1.2) without prescrib- 
ing the asymptotic behavior behind the front is largely open. The main goal 
of this paper is to provide some easily verifiable conditions to establish exis- 
tence or non-existence of multiple solutions of (1.3). To this end, we develop 



variational tools to study existence and multiplicity of supercritical fronts, i.e., 
the solutions of (1.3) with speed c > c* , where c* > is some critical prop- 
agation velocity, whose precise value will be specified in Sec. [2j One of the 
by-products of our methods is an independent proof of existence of solutions 
11 12]. However, it is worth noticing that our method 



to (1.3) obtained in 



can yield solutions of (1.3), (1.2) that are different from those of 11 12 under 



easily verifiable conditions, thus establishing multiplicity of supercritical fronts 
for certain ranges of propagation speed. At the same time, our variational char- 
acterization, together with some refined analysis of the decay of solutions ahead 
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of the front also yields global uniqueness of solutions of (1.3) and (1.2) for all 
c > c", for some c" > c* which is readily computable. Thus, the solutions ob- 



tained in 11,12 are the only positive traveling wave solutions for (1.1), (1.2 1 



invading the unstable equilibrium u = for c large enough, irrespectively of the 
asymptotic behavior at z = —oo. This is one of the main contributions of this 
paper, and we illustrate it by a simple example in which all the estimates are 
explicit. 

Our approach is mostly variational, based on the ideas developed in fT4||T8], 



relying on the fact that (1.3) is the Euler-Lagrange equation for a suitable 
functional <i>c defined on the exponentially weighted Sobolev space H^CE) (see 
Section [2] for precise definitions). However, in the application of the direct 
method of calculus of variations to get existence of the traveling wave solutions of 
interest, one runs into the difficulty that these solutions do not belong to H^{T,). 
To overcome this difficulty, we first construct, by non- variational methods, an 
auxiliary function which solves (1.3) for z sufficiently large and subtract its 
contribution from the integrand of $c, thus removing the divergence of the 
integral. We then minimize the modified functional to obtain a non-trivial 
Hl{T,) correction to the prescribed auxiliary function. As a direct consequence, 
the existence of a traveling front is established. 

We notice that the solution constructed in this way may depend on the 
choice of the auxiliary function and thus may be non-unique. Nevertheless, 
we show that this does not happen for sufficiently large values of c, in which 
case the obtained solution is the unique front solution of (1.3). The proof of 
uniqueness also has both a variational and a non-variational parts: first, using 
delicate exponentially weighted L^-estimates, we show that the difference of any 
two solutions with the same asymptotic decay at z = +oo belongs to ff^(S), 
and then we prove uniqueness by convexity of the respective functional for c 
large enough. 

The paper is organized as follows. In Section [2] we introduce hypotheses 
and notations and state main theorems. In Section |3] we construct an auxiliary 
function which represents the principle part of the traveling wave solution for 
large z. Section |4] describes the variational setting for the modified functional. 
In Section [s] we present results on existence of traveling fronts and describe 
their properties. Section |6] is devoted to the proof of uniqueness of supercritical 
fronts traveling with sufficiently large speed. We also discuss a simple one- 
dimensional example where the range of velocities for which uniqueness holds 
can be explicitly computed. Finally, in Section [7] we prove a global stability 
result for the supercritical fronts considered in Section [6] 



2 Preliminaries and main results 

In this section, we introduce the assumptions of our analysis and state our main 
results. Throughout this paper we assume f2 to be a bounded domain (connected 
open set, not necessarily simply connected) with a boundary of class C^. We 
start by listing the assumptions on the nonlinearity / which we will need in our 
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R satisfies 

/(l,2/)<0. VyeH. 



analysis: 

(HI) The function / : [0, 1] x H 

/(o,y) = o, 

(H2) For some 7 S (0, 1) 

/eC°-^([0,i] xH), /„eC°^^([0,l] xH), 

where ^ := 



(2.1) 



(2.2) 



Note that hypotheses (HI) and (H2) are, in some sense, the minimal assumptions 
needed to guarantee global existence and basic regularity of solutions of ( |1.1[ ) 
satisfying 



< u < 1. 



(2.3) 



From now on, when we speak of the solutions of either (1.1 ) or (1.3), we always 
assume that they satisfy (2.3). 

Under the above hypotheses the following functional 



is well-defined for all functions obeying (|1.2|) and (|2.3|) that lie in the exponen- 

(2.5) 



tially weighted Sobolev space H^{Y^), i.e., functions which satisfy 
/ e"^(|VMp + ii2) dx<oo, 



for a fixed c > 0. 

Our final assumption on / concerns the stability of the trivial equilibrium 
u — 0. In this paper, we consider the situation, in which this equilibrium is 
linearly unstable. This fact is expressed by the following hypothesis: 



(U) There holds 



vq := mm 

'/'lsn±=0 



In V'^ dy 



< 0. 



(2.6) 



Indeed, the constant vq defined in ( 2.6 ) is the principal eigenvalue of the problem 

Ayi;k + fuiO,y)i^k + i'ki^k=0, (2.7) 



with boundary conditions as in (1.2). In particular, for the principal eigenvalue 
we can assume ipo to satisfy 



\fpo + /«(0, 2/)?Ao + J^oV'o = 0, ■00 > 0, 



IIV'ol 



L°°(a) 



= 1. 



(2.8) 
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Note that under hypothesis (U) we have existence of a particular solution of 



(1.3) and (1.2 1 for some c* > cq, where 



(2.9) 



which plays an important role for propagation with front-like initial data 
As was shown in 



18 , there are precisely two alternatives: (i) either there exists 



a non-trivial minimizer of $(,t in H^^CS) subject to (1.2) for some > co, in 



which case c* = c^; or (ii) there are no non-trivial minimizers of <i>c in H^iYi) 
subject to ( 1.2 ) for all c > cq, in which case c* = cq and there is a minimal wave 
with speed cq. For our purposes here, it is convenient to introduce the following 
equivalent characterization of c*: 

c* inf |c> 



inf $cM = 

0<M<1 



(2.10) 



where the infimum of <i>c is taken over all functions in HI{Yj) satisfying (1.2) 
The equivalence of the two definitions of c* is shown by Proposition 5.1 
For c > c* any positive solution of (1.3) and (1.2) is expected to satisfy 



(2.11) 



^I(y,z)~aVo(y)e-^-^^^''"'^ A±(c,i/) 

as z — -|-oo, with some a > and with -00 > being the minimizer of the 
Rayleigh quotient in (2.6) normalized as in (2.8). At least formally, the asymp- 
totic behavior of positive solutions of ( 1.3 ) and ( |1.2[ ) can be obtained by lineariz- 
ing (1.3) around w = 0, leading to the formula in (2.11) with the exponential 



decay rate determined by either A+(c, I'd) or A_(c, fo). However, the case cor- 
responding to A+ (0,1^0) is impossible, since then the solution would belong to 
i?^(S), i.e., it would be a variational traveling wave (for the definition and fur- 
ther discussion, see [T6}|l8) ). But by [T7[ Proposition 3.5] this contradicts the 
assumption c > c*. In other words, for c > c* the traveling wave solutions 
acquire fat tails. 

It is also well-known that solutions of ( |1.3[ ) become z-independent as z — > 
—00. Thus, as z — > —00 the function ?2(-,z) is expected to approach a critical 
point of the energy functional 



E[v] 



-\Vyv\' + V{v,y)]dy, 



(2.12) 



which is defined for all v E H^{il) with values in [0,1] satisfying Dirichlet 
boundary conditions on 9ri±|^ By hypothesis (112) the critical points of E 
satisfy the equation 



+ /(w,y) = 0, v\on^ 



0, i/-Vw| 



9^0 



0. 



(2.13) 



^In the variational context throughout this paper, the Dirichlet boundary conditions are 
understood in the usual sense of considering the completion of the family of smooth func- 
tions whose support vanishes at the Dirichlet portion of the boundary with respect to the 
corresponding Sobolev norm. 
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Similarly to the case of z — > +co discussed in the previous paragraph, the 
convergence of the traveling wave solution to v = liuiz^-oo z) satisfying 
(2.131 can be, at least formally, obtained from the analysis of (1.3) linearized 
around v. Introduce 

. U\'^y^\'-fu{v,y)i'')dy 
I'o = mm — — = , (2-14) 

whose principal eigenvalue satisfies 

Ayi/io + /«(0,2/)Vio + i>oV'o = 0, Vio>0, ||V'o||L~(f2) = 1- (2.15) 

Then in the generic case f'o > Oi when u is a strict local minimizer of E, 
one expects the solution to approach v from below when z — > — cx3 as 

?2(y,z) ~ t>(y) + ael^-^'^''^''^l^-0o(y) for some 5 < 0, (2.16) 

where we noted that the rate of convergence is determined by the only negative 
value of A±(c, vq) [Tsj. 

We now turn to our result on existence of solutions to (1.3) and (1.2) that 
admit a particular variational characterization. These are the traveling wave 
solutions which are minimizers of the functional $™ defined as 

<I)-M=^e-Q|V^x|2 + F(u,y)-i|Vw|2-T/(i«,y)^ dx, (2.17) 

for some suitably chosen w : E — K. Specifically, the function w e C^(E) n 
C"'^(S) should solve (1.3) and (1.2 1 far ahead and far behind the front. As we 
show in Sec. |4] the functional is then well defined for all u satisfying (2.3 1 

and u-we Hl{T,). 

Our existence and variational characterization result is given by the following 
theorem. 
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Theorem 1. Assume hypotheses (HI), (H2) and (U) are true. Then, for every 
c> c* andaGR there exists u S C'^{T.)'r^C^(Y.)f^W^'°°{T.) such thatO <u<l 
and: 



(i) u solves (1.3) and (1.2), and there exists w G C (S) H C (S) with values 



in [0, 1] solving (1.3) and (1.2) in fix (— oo, — Zq) U x (zq, +(X)) for some 
Zq > 0, such that the function u minimizes <i>™ [u] over all u satisfying 
([t3|, such that u-w e H^{^). 



(ii) u{y,z) — a^ Q(y)e •^-(^^'^0)2 q^^ Az^^ some A > A_(c, I'o) and '00 
solving (2.8), uniformly in C^ijl x [i?, +00)), as R ^ +00. 



(Hi) u is a strictly monotonically decreasing function of z, and we have u{-, z) — >■ 
V uniformly in C^{Q.), where Q < v < \ is a critical point of E, and 
E[v] < 0. 

(iv) I'o > 0, moreover, if > 0, then u{y,z) = v{y) + d'ipo{y)e~^-^''''^"^^ + 
0(e^'^^), with some a < and A < A_(c, ;>o), uniformly in C^{il x 
(—00, R]), as R ^ —00. 



(v) If w e C^(S) n C^(S) with values in [0, 1] solves (1.3) and 
(—00, —Zq) U r2 X (zq, +00) for some zq > and w ~~ w ^ 
also minimizes 



1.2) in CI X 



(E) , then u 



Let us point out that existence results for supercritical fronts in cylinders 
were first obtained in the classical work of Berestycki and Nirenberg under a 
few extra assumptions [IT] (see also [l2{[T3] ). Specifically, they assumed (in our 
notation and setting) that if vi is the smallest positive critical point of E (in 
the sense that if < u < 1 is a critical point of E, then v > vi) and if vi is 
stable, i.e., if i^o > with u = t;i in (2.14), then there exists a unique (up to 
translations) traveling wave solution approaching zero from above as z — ?■ +00 
and vi from below as z — > —00. Note that existence of vi and the property 
I'o > for V = vi follows immediately from our hypotheses. 

It is easy to see that the existence and monotonicity result for fronts con- 
necting zero and vi may also be obtained from a variant of our Theorem [l] 
in which a minimizer of $!f is constructed in the class of functions bounded 



above by vi. More precisely, let us define the value of cj by (2.10), in which the 
infimum is taken over all u € Hl{T,) satisfying (1.2) and such that u<vi: 



= inf <^ c> 



inf $cM = 

0<u<Di 



(2.18) 



where, again, the infimum is taken over all u G iJ^(E) satisfying (1.2 1. Then 



under hypotheses (HI), (H2), and (U), for every c > > Cq and every a € E 
there exists u which satisfies the conclusions of Theorem [T] with v = vi, where 
vi is the smallest positive solution of ( |2.13 1. In particular, it can be seen from 
Theorem [1] that u is the minimizer of $Jf for all w with values in [0, 1], such 
that w-ue Hci^)- 
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It is important to point out that using their sUding domains method, Beresty- 



cki and Nirenberg also proved uniqueness of solutions of (|l.3| and (|1-2|) con- 
On 



necting zero and a non-degenerate vi for all c > c| (sec also 12 13 



19 



the other hand, a priori one cannot expect uniqueness in the more general class 
of solutions satisfying ( |2.3[ ), even in one dimension (see the counter-example in 
Sec. [g]). a natural question, then, is whether for a given c > c* the minimizer 
given by Theorem [T] is the solution obtained by Berestycki and Nirenberg. We 
can give a negative answer to this question in the case when c* > for a certain 
range of c > c*. 

Indeed, suppose c* > c^. Then, since > Cq, there exists a non-trivial min- 
imizer u* of $c* among all < u < 1. Furthermore, since u* is not a minimizer 
of O'^sr all u < vi, we have lhnz-f-oo'u*{-, z) > vi. In particular, this 

means that for 



inf <^ c> : inf [u] = 



Vi<U< 1 



(2.19) 



we have c|,j^ > c*. If not, for ui = min(u*,wi) and U2 = ma.x{u* ,vi) we would 
obtain = 'i'c* [^i] + '^'c^[w2] > *&c*[^i]' implying that ui < vi is a non- 

trivial minimizer. However, this contradicts the fact that there are no non-trivial 
minimizers of $c* that lie below vi . 

In fact c|, > c*, since the infimum in (2.19) is attained, with the value of 



cl^ being the speed of the non-trivial minimizer of ^""^ over all u, such that 

u — vi E H^CS) and vi < u < 1. Indeed, there exists a sequence c„ — > cl_^ from 
below and a sequence of functions such that [un] < 0. Therefore, 

in view of the fact that i>o > toi v = vi, by |18j Theorem 3.3] applied to ^'^^ 
(see also the proof of 20 Lemma 3.5]), there exists a non-trivial minimizer of 
4>"i for some c > c ] , a nd by the definition of cj,^ and the arguments in the 
proof of Proposition 



5.1 



we have c = c. 



Since by the argument in the proof of 20 Lemma 3.5] for every < m < fi, 
such that u — w€ if^(I]), with w as in Theorem[l] and u{-, z) vi in C^{fl) as 
z — — oo, the value of $^[m] can be lowered by extending u above vi for every 
c < cl , we have the following consequence of Theorem [ll 



Corollary 2.1. Let c* > c*^ 

solution u in Theorem^ has v > vi. 



Then cl^ 



> c* , and for every c* < c < cl the 



In other words, if c* > c^, for every c as in Corollary 2.1 



exist at least two solutions of (1.3) and (1.2) satisfying (2.111 



and every a G M there 



Remark 2.2. Note that the argument leading to Corollary 2.1 also works under 
the assumption that > c* when c* = cl = cq. 

Let us now go back to the question under which conditions the solutions 
of Berestycki and Nirenberg, which connect zero with vi for every c > are, 
in fact, the unique (up to translations) solutions of (1.3 1 and (1.2). In one 
dimension it is possible to use phase plane arguments to show that for large 
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enough values of c the only solution of ( 1.3 ) connects zero to vi (e.g., under non- 
degeneracy of wi). What about the higher-dimensional case? In the following we 
show that the solution connecting zero to vi is the only traveling wave solution 
for large enough speeds, up to translations. 
We start by introducing 



vq '■= nim 



In V'^ dy 



liy) ■= max A(s,2/). 
se[o,i] 



In analogy with (2.9 1, we can then define the quantity 



(2.20) 



(2.21) 



which is easily seen to satisfy cq < c" < oo. It turns out that gives a lower 



bound for the values of c, for which the solution of (1.3) and (1.2) is unique. 



up to translations. Our global uniqueness result, which relies on the variational 
characterization of the obtained solutions is contained in the following theorem. 



Theorem 2. Let c > , where > Cq is defined in (2.21), and let the assump- 
tions of Theorem [7] be satisfied. Then > c* , and for every c > there exists 
a unique solution of ( |1.3[ ) and (1.2) (up to translations). 

Furthermore, the unique solution of Theorem [2] is globally stable with respect 
to perturbations with sufficiently fast exponential decay. 

Theorem 3. Let u he the traveling wave solution with speed c > from The- 
orem^ and let < Uq < 1 satisfy Uq — w G L^, (S), for some c' > 2A_(c, Vq). 
Then, if u is the solution of (1.1^ and {1.2) with u{x,0) — uq(x), we have 

(2.22) 



||u(?/,z + ct,t) - u(y,z)||i2^(s) 



for some cr > 0, C > 0. 

We note that the result of Theorem [3] also holds in other norms, e.g. with 
respect to uniform convergence on compacts (see 20 for details). 

Note that as a direct consequence of Theorems [l}|3] we have a complete 
characterization of the supercritical fronts in problems involving nonlinearities 
of the kind that was originally considered in [T0II2T] : 



Corollary 2.3. Under hypotheses (HI), (H2) and (U), assume that, in addi- 
tion, fu{u,-) < /m(0, •) for all u £ [0,1]. Then for each c > Cq there exists a 
unique (up to translations) solution of ( |1.3[ ) and \1.2\ , which is globally sta- 
ble in the sense of Theorem^ with respect to perturbations in i^, (E), with any 
d > 2\_{c,va)- 

This corollary is obtained by noting that under its assumption on the nonlin- 
earity we simply have — uq. Therefore, c" = cq and, in particular, c* — cq as 
well. In fact, in view of the conclusion of Theorem [ijii) the following conver- 
gence result holds under a single assumption on the rate of exponential decay 
of the initial data. 
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Corollary 2.4. Under the assumptions of Corollary \2.3[ let c > cq and assume 
that uoiy, z) = a-0o(y)e"'*'-'^'''''°^^ + 0(e~'^''), for some a e M and X > A_(c, i/q). 
Then the conclusion of Theorem^ holds. 

In the rest of the paper (HI), (H2), (U) are always assumed to be satisfied. 



3 The front far ahead 



In this section we start the investigation of supercritical traveling wave solutions. 
As was mentioned in the introduction, the first important step in the study of 



traveling fronts is the construction of an auxiliary function which solves (1.31 



and (1.2) for z large positive and has the asymptotic behavior given by (2.11) 



Specifically we will construct a function w with the following properties: 
(Wl) w e C2(E) n Ci(S) and solves 

Aw + cw^ + f{w,y) = 0, 



(3.1) 



in S for z > 3, with the same boundary conditions as in (1.2) 
(W2) There exists 5 > and a > sufficiently small, such that 

w{y, z) = aV'o(y)e~-^-^"''"'^^(l + 0(e"''")), z ^ +oo 



(3.2) 



where A_(c, z^g) £^nd ■00 are defined in (2.11) and (2.8). 
(W3) We have 

< w{v, z) < 2aMy) < 1 V(y, z) e S, 
and u> = for all z < 2. 



(3.3) 



In order to establish existence of a function w satisfying (Wl) - (W3), we 
only need to prove existence of a solution of (3.1 ) satisfying (W2) and (W3) for 
z > 2. Then we obtain the desired function by multiplying this solution with a 
cutoff function 77 G C°°(]R), such that 77(2) = 1 for z > 3 and ri{z) = for z < 2. 

By translational symmetry in the z-direction, any translate of solution of 



(3.1) is once again a solution for sufficiently large z. Translations of w are 



equivalent to adjusting the value of the constant a in (3.2). Therefore, fixing 



the constant a in (3.2) is equivalent to fixing translations of the solution of 
(3.1). Thus, the choice of a particular value of a is inconsequential. Indeed, 



considering solutions with a fixed and z large enough positive is equivalent to 
considering the problem for small enough a > and z > 0. Conceptually, this 
choice of w resembles selecting one element of the stable manifold corresponding 
to the fixed point m = in the one-dimcnsional setting (in the sense of ordinary 
differential equations). 

The following proposition establishes existence of a function w with all the 
desired properties. 
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Proposition 3.1. For every c > cq and every a > sufficiently small, there 
exists a function w with properties (W1)-(W3). 

The proof of this proposition is based on two lemmas. We first explicitly 
construct sub and supersolutions (see definitions (i) and (ii) below) for the 
problem (3.1) on the right half cylinder (Lemma |3.2[ ). Next we show that if 
one restrict the problem (3.1) on bounded sections of the cylinder, then this 



problem has a solution squeezed between sub and supersolutions (Lemma 3.3) 
We then complete the proof of Proposition |3 . 1 1 by passing to the limit. 



In order to proceed, we introduce the notions of sub and supersolutions: 



(i) We call w a supersolution of the problem (3.1), if w G C (il x (0, +oo)) 



satisfies the boundary conditions in (1.2) and 



Aw 



cw. 



/(W,j/)<0, z>0. 



with 



w(0,j/) = 2a-0o(2/), y&^- 



(3.4) 



(3.5) 



(ii) We call w a subsolution for the problem (3.1), if w G C^{Vl x (0, +oo)) 
satisfies the boundary conditions in ( |1.2[ ) and 



ith 



Aw + cw^ + /(ui,?/) > 0, z > 0, 
«Z(0,2;) = 0, yen. 



(3.6) 
(3.7) 



Lemma 3.2. There exist w and w, the sub- and super-solutions of (3.1 ), which 
satisfy (W2) and (W3), for every a > small enough. 



Proof. The proof follows by explicit construction of w_ and w verifying (1.2), 



( |3.4[ ) - (3Jj), (W2), and_(W3). By assumption (H2), for ||w||l~(s) < 2a < 1 
there exist constants fc, fc > 0, such that 



fu{Q,y)w - k\w\'+'< < f{w,y) < fu{0,y)w + k\w\'+\ 



Therefore, the function w satisfying (1.2), ( |3.5[ ) and 

7V(w) Aw + cw^ + /„(0, y)w + kw^+'' < 0, 



in 51 X (0, +oo) automatically satisfies (3.4), and thus is a supersolution. 
Let us show that the function 



w{z,y)^aMy)e'^-'{^ + e''l 



-,-Sz 



(3.8) 
(3.9) 

(3.10) 



is a supersolution for some < (5 ^ 1 (throughout the proof we use the short- 
hand A_ for A_(c, i^o)). Indeed, substituting this function into (3.9|, we obtain 



U{w) = a{Ay^Po + i^--cX-+fuiO,yWo)e-^-^l 
-a6{c - 2A_ - S)tPoe^^^-^^^" 



(3.11) 



l+7,/,l+7^-(l+7)A-Z 



(1 
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The first line in (3.11 1 is identically zero in view of the definition of A_ and V'Oj 
see and ( |2.11[ ). Therefore, 

-a<5(c-2A_ -(5)V'oe~(^-+*^^ 



S{c- 2X^-6) 



(3.12) 



Since also c > 2A_ by (2.11 1, the factor multiplying the expression in the curly 
brackets in (3.121 is negative for S sufficiently small. Let us show that the 
expression in the curly brackets can be made positive by a suitable choice of a 
and S. Indeed, since the expression in the curly bracket is an increasing function 
of z for i5 small enough, we only need to verify its positivity at z = 0. Then, 
the result easily follows by choosing S sufficiently small first, and then choosing 
a such that a'' 5. 

We now turn to sub-solutions. By the same argument as above, it is enough 
to show that 



M{w) Aw + cw^ + fu{0,y)w-kw^+^ > 



in 51 X (0, +oo) for some w>0. Similarly to (3.101, we choose 



5z\ 



(3.13) 



(3.14) 



Substituting this into (3.13) and using the same arguments as for the superso- 
lution, we obtain 



Miw) = aS{c^2X^ -6)^06-'-^-+^^'' 



1 - 



S{c- 2X^-5) 



(3.15) 



which is positive for a similar choice of S and a. 

Finally, we note that, as a consequence of our construction, the functions w 
and w satisfy (3.2) and (3.3). □ 



We next establish existence of a classical solution of (3.1 ) on finite sections 
of the cylinder away from the boundaries. 

Lemma 3.3. LetR>l. Then there exists wr e {fix {1, R))r\C^{Tlx[l, R]) 
solving (3.1) in Q x [1, R], and satisfying w < wr < W, where w and w are as 
in Lemrna \3.S\ 



Proof. First, the proof of existence of a weak solution of (3.1 ) in S/? = fix (0, R+ 
1), sandwiched between the sub- and the super-solutions of Lemma 3.2 follows by 
standard monotone iteration argument (see, e.g., 22 Sec. 9.3].) Furthermore, 
by standard elliptic regularity, this solution is regular in the interior of [23| . 
Now, we apply regularity theory to the solution in E/; = il x (1, i?) CC E/j 



23 Theorem 9.13] to obtain regularity up to the boundary on the Dirichlet 

< C"||wfl||LP(SH): (3-16) 



portions dT,± of the cylinder: 

\\wR\\■^^2,p^^^) < C{\\wr\\lp(^^^) 
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for some C,C' > and any p < oo, where we took into account hypothesis (H2) 
on /. Hence, by Sobolev imbedding wu is in C^'"(Sfl), for any a e (0, 1). By 
construction, w < wr < w in the closure of S^. □ 



Proof of Proposition 3.1. We now finish the proof of Proposition 3.1 by passing 
to the limit i? — > oo in Lemma 3.3 Observe that by the bounds of Lemma 3.3 

we have 



and the exponential decay of w; and w of Lemma 3.2 



\WR\\LP{n y.{l,R )) 



< 



C, where C > is independent of R. Therefore, by the estimate of ( 3.16| ) we 



have wr 



Wy^^{VL X (l,oo)) on a sequence of 



In view of 



arbitrariness of p < oo, the convergence of wr to w is al so in C^{Q, x [2, i?— 1]), 
for any i? > 3. Hence, the limit w is a weak solution of (3.1 ) in x [2, oo), and 
by standard elliptic regularity is also a classical solution. The proof is completed 
by multiplying the obtained solution by a cutoff function. □ 

Remark 3.4. Our proof guarantees existence of infinitely many functions sat- 
isfying (Wl) - (W3), since existence of a solution for 
fixed boundary data at z — satisfying (3.3). 



3. 1 ) also holds for each 



4 Variational setting 



As was already noted in the Introduction, certain solutions of may be 

viewed as critical points of the functional $c[w] introduced in ( |2.4[ ). This, how- 
ever, poses some restrictions on the type of traveling wave solutions (which are 
the so-called variational traveling waves 16 18]) captured by this approach. In- 
deed, in view of the presence of the exponential weight in the functional the 
set of admissible functions must be characterized by the natural exponential 
decay at z = +oo. At the same time, it is known (see, e.g. 11 13]) that when 
u = is an unstable equilibrium of (1.1 1, generically there are traveling wave 



solutions which fail to have this natural exponential decay. In this situation 
the use of a variational approach to the construction of traveling wave solutions 
seems problematic due the divergence of the respective integrals, even though 
formally (1.3) is still the Euler-Lagrange equation for We propose to elim- 



inate this divergence by subtracting the divergent part of the integrand in the 
definition of Namely, we start by introducing the modified functional in 

(W3) of Sec. H In 
its value 



(2.171, where u; is a function satisfying properties (Wl) 
other words, we subtract from the integrand in the definition of <I>c 



obtained by substituting an exact solution of ( 1.3 1 for z > 3 (property (Wl) of 
w) and characterized by a prescribed exponential decay at z = -|-oo (property 
(W2) of w). Note that existence of w with these properties was proved in Sec. 
[3] by non- variational techniques. 

We next identify the suitable admissible class of functions u on which the 
new functional [w] makes sense. Let u — w + h, for some unknown function 
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/i : S — > M. With this definition, we now have = where 

'1 



E 



*rW= / e''[;^\Vh\^ + Viw + h,y)-V{w,y)-V'iw,y)hjdx 

Aw + cw^ -V'{w,y)]hdx, (4.1) 



e 

E 



e 



where here and throughout the rest of this section we use the notation V'{u, y) := 
dV{u,y)/du and V"{u,y) := d'^V{u,y) / dv? . Observe that as long as u{x) S 
[0, 1] for all a; G S, we can use Taylor formula to rewrite vE'™[/i] in the following 
form 

= ^ Q iV/ip + \v"[u, y)h^^ dx 

Aw + cWz-V'{w,y)^hdydz, (4.2) 

for some u sandwiched between w and u. Here we took into account hypothesis 
(H2) and that, by property (Wl) of w, the integrand in the last term in (4.11 
is identically zero for z > 3 and z < 2. Now, since V' ,V" ,w,Wz are uniformly 
bounded and Aw £ L'^{n x (2, 3)), the functional vjfjf' is well-defined within the 
admissible class , defined as follows: 

Definition 4.1. A function ft- : E — !■ M is said to belong to the admissible class 
A^, if he Hl{T) and < h{x) + w{x) < 1 for all a; € E. 

Note that the assumption that h G implies that u = w + h has values in 
the unit interval is not too restrictive, in view of hypothesis (HI), since u = 
and u = 1 are sub- and super-solutions of (1.3 1 and, hence, are the natural 



barriers 24 



In order to apply the direct method of calculus of variation, we need to 
establish weak sequential lower-semicontinuity and coercivity of the functional 
\E'™ in the considered function class. Lower-semicontinuity in the considered 
class of problems was studied in 17 , and we have 

Proposition 4.2. For all c > cq, the functional is sequentially lower- 
semicontinuous in the weak topology of H^(Y^) within A^ . 



Proof. The proof follows exactly as in 17 Proposition 5.5]. □ 

Thus, the main difficulty in applying a variational approach to supercritical 
fronts has to do with coercivity of ^Jf. Below we show that coercivity indeed 
holds for all c > c* . 

Proposition 4.3. The functional is coercive in A^ for every c > c* . 
Proof. We carry out the proof in 4 steps. 
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step 1 First of all, since h, V , w, Wz are uniformly bounded and since Aw G 
L^{rt X (2,3)), by a straightforward application of Cauchy-Schwarz inequality 
we have 



2 Jn 



Aw + cwz — V [w, y) I h dydz 



(4.3) 



for some Ki > 0. On the other hand, using the formula in the proof of 17 



Proposition 6.9] and the fact that by the definition of c* we have $c' [p] > for 
all c' G (c*, c) and ah p G Hl,{T,), it holds that 



Therefore, we have the following estimate: 



2 /2 

— c 
2c2 



(4.4) 



+ £ e^^ (^F(u; + h,y)- V'{w, y)h - ^(u;, y) - V{h, y)^ dx 



(4.5) 



> 



2 /2 

C - C 



2c2 



+ ^ e^^ (u; + /i, y) - V'{w, y)h - y (w, y) - V{h, y)^ dx. 

Step 2 By Taylor formula, there exist w and ft, such that — < \h\ and 
\h\ < \h\ such that 



V{w + h,y)- V'{w, y)h - V{w, y) ~ V{h, y) 



]^{v"{w,y)-V"{h,y) 



On the other hand, hypothesis (H2) implies that 

\V"{w,y)~V"{h,y)\<C^w-h\\ 



(4.6) 



(4.7) 



for some Ci > 0. Therefore, for all \h\ < a, where a is defined in property (W3) 
of w, we obtain 



\V(w + h,y)- V\w, y)h - V{w, y) - V{h, y)\ < C^a^^h^ 

for some C2 > 0. 

Now, by boundedness of V' , V" , when \h\ > a we also have 



(4.8) 



\V{w + h,y)~r{w,y)h~V{w,y)-Vih,y)\<C3W, (4.9) 
for some C3 > 0. 
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step 3 Using (4.6)-(4.8) and the Poincare inequality 17 Lemma 5.1], one 
can write 



{|/i|<a} 



e^^l V(w + h,y)- V'{w, y)h - V{w, y) - V{h, y)\dx 



(4.10) 



Similarly, using (4.9) and property (W3) of w we have 

/ e'=^|F(u; + h,y)- V'{w,y)h - V{w,y) - V{h,y)\dx 

J{\h\>a} 
R 



e'^^'wdydz + C3 



-00 J {\h{- ,z)\>a} 



R "'{|h(-,2)|>a} 



e'^^w dydz 



^cR 



R J{\hi-,z)\>a} 



e'^^'wdydz, (4.11) 



for any R > 0, where denotes the (n — l)-dimensional Lebesgue measure of 
ri. Then, since by Poincare inequality |17l Lemma 2.1] we have 



e'<hldydC> / h^y,z)dy 



> 



{\h{-,z)\>a} 



h^{y, z) dy > 



dy. 



{\h(-,z)\>a} 



(4.12) 



Property (W2) implies that w{y,z) < 2ae for sufficiently large z. Using 
this observation and inequality (4.12), we obtain 

( p(c-A-)2 



R J{\h{-,z)\>a} 



e'^^w{y, z)dydz < 2a 



dy dz 



{\h(-,z)\>a} 



< 



ac 



R 



-\^z 



f f e'^'^hldydndz 
Jz Jn J 



2 

< — 
ac 



-^-''dz 



ac\- 



R JQ. 



e'^^'hl dydz < 



^czu2 

R Jn 

A__R / czi2 



acX- 



e'^'h'l dydz 
e^^hldx.iA.U) 



It then follows 



{W>a} 



e^^ \V{w + h,y)- V'{w, y)h - V{w, y) - V{h, y)\dx 



<Ciae''" + C5a-^e-^-" [ e^'hldx, (4.14) 



for some 6*4,(75 > 0. 
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In view of the above, for a given a one can choose R sufficiently large so that 
e""" \V{w + h,y)- V'{w, y)h - V{w, y) - V{h, y)\dx 



m>a} 



where K2 = K2{a) and Cg > is a constant independent of a. 



(4.15) 



Step 4 Finally, combining the estimates of Steps 1 and 3 and choosing a small 
enough, we obtain 



> -Kj, + Cric -c* -e) [ e'^'hl dx, 



(4.16) 



where i^^3, C7 > 0, for any given e > 0. Therefore, choosing e sufficiently small 



and also making use of the Poincare inequality 17 Lemma 2.1], we arrive at 

-''^h^dx<K, (4.17) 



e^'hldx < K, 



for some K > Q on a. sub- level set of ^t™. This, in turn, implies that 
/ e""|Vj,/ipda; < 2*;f[/i] + j e''''\V'\u,y)\h'^dx < K\ 



(4.18) 



for some K' > 0, on that sub-level set as well, and hence the functional is 
coercive for all c > c*. □ 



5 Existence and properties of minimizers of $ 



In this section we establish existence of traveling wave solutions of (1.1) with 



exponential decay governed by (2.111, which are global minimizers of the re- 



spective functionals. Let us point out that in general this approach may not 



give all such solutions of (1.3) and (1.2 1. Nevertheless, as we will show later in 
Section [6] these are in fact the only such solutions for large enough values of c. 
Also, for certain classes of nonlinearities, the same result holds for all traveling 



wave solutions with c > cq (see Corollary 2.3 1 



We begin with a basic characterization of c* in (2.10) in terms of existence 
of minimizers of the functional $c for c > co- 



Proposition 5.1. Let c* be defined in ( |2.10 |. Then c* > cq, there is no non- 
trivial minimizer of $c over all u G HJ(T^ satisfying (1.2) for any c > c* , 
and if c* > cq, there exists a non-trivial minimizer of ^c* over all u G iJ^*(E) 



satisfying (1.2) 
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Proof. Suppose first that there exists a non-trivial minimizer of $c for some 
c = c' > Co. Then it is clear that c* = . Indeed, by (l?! Proposition 3.2] and 
the arguments of \17\ P roposition 6.9] we have inf $c = for every c> . On 
the other hand, by [18[ Equation (5.4)] we have inf $c < for all < c < c^. 

Alternatively, if there is no non-trivial minimizer for any c > cq, then 
inf$c = for all c > cq. Indeed, if not then there exists a trial function 
Uc € ffc(S) with values in [0,1], such that <I>c[uc] < 0, contradicting non- 
existence of minimizers guaranteed by 18 Theorem 3.3]. On the other hand. 



arguing as in [l8j Theorem 4.2], it is not difficult to see that the trial function 
in 18 Equation (4.18)] makes $(, negative for any < c < Cq. 



□ 



The proof of existence of traveling wave solutions in the statement of Theo- 
rem [l] follows from a sequence of propositions below. 

Proposition 5.2. Let c > c* and let w satisfy (Wl) - (W3). Then there exists 
u e C2(I]) n Ci(S) n Vri'°°(S), with O < w < l, such that 



(i) u solves (1.3) with 



(a) u minimizes '^^[u] over all u taking values in [0, 1] and such that u — w G 



Proof. By coercivity of established in Proposition 4.3 if {h„) is a minimizing 
sequence from , then (/i„) is uniformly bounded in iJ^(E). Hence there exists 
a subsequence, still labeled (/?,„) such that hn ^ h in H^{I]). Moreover, h is 
obviously still in A^' . Therefore, by low er-s emicontinuity of ^'^ in the weak 
topology of i?^(S) proved in Proposition 
in A^. Now, since [u] = '^f[u — w], the function u — 



the limit /i is a minimizer of 

h is also a minimizer 



of <I>Jf' over all u, such that < m < 1 and u — w G H^(I]). This proves part (ii) 
of the statement. 

To prove part (i) , consider the first variation of the functional $^ evaluated 
on u with respect to a smooth function ip with compact support, vanishing at 
dT,±. Since u is a minimizer and the barriers u ~ and u — I are sub- and 



supersolutions of (1.3), we have 24 



e"'' {Vu • Vip - f{u, y)(p} dx. 



(5.1) 



23 M is a classical solution of (1.3), and u £ 



= 5$^(u,(^) = 

By elliptic regularity theory 
M^i'°°(S). 

Moreover, we also have u ^ 0, since u — would imply h — u ~ w ^ —w G 
-ff^(E), contradicting our assumptions on w. By strong maximum principle, we 
then have < u < 1, in view of the fact that u = and u = 1 are the sub- and 
super-solutions of (1.3 1, respectively 18 23 . □ 



Remark 5.3. We note that the proof of Proposition \ 5.S\ does not rely on the 
precise details of the assumptions (W1)-(W3) on the function w. The main 
ingredient of the proof is that w solves (|1.3[) and ( 1.2 ) and does not lie in H^{T,). 
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Therefore, with a suitable construction of w our approach may also be applied 
to supercritical traveling waves in the case vq — 0, yielding existence of traveling 
wave solutions with non- exponential decay. 

It is not difficult to see that the minimizer constructed in Proposition [5?2] has 



the decay governed by (2.11). We note that exponential decay of solutions of 



( 1.3 1 in cylindrical domains has been extensively studied under various assump- 



tions on the nonlinearity, etc. (see, for instance, 11 13]). In particular, with 
minor modifications the arguments of 13 could be applied in our slightly more 



general setting. For the sake of completeness, we present the proof under our 
set of assumptions, using the techniques of ,18j and giving a slightly stronger 
decay estimates. We also note that our result establishes the precise rate of 



exponential decay given by (2.11) for all solutions of (1.3) with c > c*. 



Proposition 5.4. Let c> c* and let u be a solution of (J. j| ) and (1.2). Th 



u{y,z) = ae-^-(='''o)>o(y)(l + 0(e-*'^)), (5.2) 
uAy,z) = -aA_(c,i.o)e-^-("''^°)"^o(2;)(l + 0(e-*^)), (5.3) 

uniformly in Q for some a € M and S > 0. 

Proof. We reason as in the proof of 18, Theorem 3.3(iii)], recalling that c > cq 
by Proposition 5.1 By a Fourier expansion of (|1.3|) in terms of ipk and i^k from 



(2.7), after some algebra we obtain that ak{z) — ^^u{y,z)'ipk{v)dy G C^(M) 
satisfy 



COu 



Vk^k = 5fc, 



where 



9k{z)^ / {fu{0,y) - fu{u{y,z),y))u{y,z)%pk{y)dy, 
Jn 



(5.4) 



(5.5) 



for some < u < u. Using variation of parameters, for any i? e M we can write 



the solution of (5.4) in the form 



1 

Vc^+4:l^k JR 
1 

+ 4i/fc 



„A+(c,jyfc)(^-z) 



(5.6) 



where X± are defined in (2.11 ). Specifically, since ipo > and m > 0, for fc = 
we have oq > and 



\goiz)\ < C\\u{-,z)\\l^f^^^ao{z) 



(5.7) 



for some C > 0. 
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Since z) — >■ uniformly as 2: — ^ +00, by the definition of 00(2:) and (5.71 
we have 



00(2:) — ^ and 



9q{z) 
ao{z) 







-00. 



(5.; 



Therefore, by (5.6) for every 60 > there exists Rq £ M, such that for every 



6 G (0, 60) and every Rq > Rq we have 

ao(i?o) < a([(i?)e^+ + ao(^)e^"'''""^' + ^ max aoiz), (5.9) 

z Ro<z<R 



where R ^ Rq + b. In fact, the maximum in the last term of (5.9) is attained 



at 2: = ii!o due to monotonic decrease of ao{z) for all z > Rq, provided that Rq 
is large enough. Indeed, in view of (5.8) we can always choose ^0 sufficiently 
large, so that max^>;jjj ao{z) = ao(-Ro)- Then, if 09(2) is not monotonically 
decreasing, it attains a local minimum at some zq G (Rf),+oo). However, by 



(5.4) with k = and (5.8) we have 



ao(^o) < h'^oaoizo) < 0, 



(5.10) 



contradicting minimality of ao{z) at zq. 

Recalling that b = R — Rq, the monotonicity of ao{z) for z > ^0 and the 
estimate in (5.9) imply that for every e > there exists 5i > independent of 
Rq, such that 



a+{R)- 



rao(i?) Vi?e [i?o + 6i,i?o + 26i], 



(5.11) 



where the last line follows from A+(c, i^o) > ^-{c,i^q) > 0, choosing bo large 
enough. Therefore, by (5.6) and, in particular, the fact that ao(-R) — cioiR) + 
aQ{R) we have 



dao 
dz 



= -A_(c, t/o)ao (i?) - A+(c, i^o)aj(^) 
1 



< 



-A_(c,z/o)ao(i?) < Vi?e [i?o + 6i,i?o + 2&i],(5.12) 



provided that e is small enough. Now, in view of arbitrariness of -Rq ^ Rq-, 
we conclude that the inequality (5.12) in fact holds for all z > Ri := Rq + 61, 
and so 00(2) < ao(i?i)e-5-^-('=^'"o)(^"-"i) for aU z > Ri. Hence, by 18, Eq. 
(4.24)] we have u{-,z) < Ce"''^, for some C > and ^ > 0. From now on, 
we can argue exactly as in the proof flS) Theorem 3.3(iii)] to show that either 



u{y,z) 



V'o(2/)(l + e ^'') or u{y, z) = ae 



i/^o(2/)(l + e-*") 



for some a > and S > 0. Here we took into account that the slightly stronger 

follows from the 
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multiplicative decay estimate than the one obtained in 
convergence on slices of S for the difference between u and the leading order 



term in 18 Eq. (3.38)] and the Hopf lemma applied to iJjq. 
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To finish the proof of (5.2 1, we need to exclude the possibility that the 
exponential decay rate is governed by A+(c, uq). Indeed, if the latter were true, 
then by the above estimates together with similar estimates for |Vm| following 
from (18 



Eq. (3.38)], we find that u G 7J^(E). Therefore, u satisfies hypothesis 
(H3) of 18], implying that there exists a non-trivial minimizer of ioi some 
c' > c. Then, by the arguments in the proof of Proposition 5.1 this would imply 



that c* > c', thus leading to a contradiction. 

Finally, the estimate for Uz in (5.3) is obtained by differentiating (1.3) with 
respect to z and using the same arguments as for (5.2). □ 



Thus, according to Propositions |5.2| and |5.4| there exist bounded solutions 
3) and (1.2) with the behavior at z = +qo governed by the slowest expo- 
nentially decaying positive solution of the respective linearized equation. Now 
we will establish a number of additional properties of these solutions, based on 
the minimality of ^'^ evaluated on u. In particular, we will show that these 
solutions have the form of advancing fronts. 



Proposition 5.5. Let u be as in Proposition \5.S\ Then 

(i) u{y,z) is a strictly monotonically decreasing function of z. 

(ii) u{-,z) ~> V in C^{il) as z ^ — oo, where < v < 1 is a critical point of 
E, and E[v\ < 0. 



Proof. We prove monotonicity by using a monotone decreasing rearrangement 
of u{y, •), following an idea of Heinze 15 . Introducing a new variable C = e'^^, 
we can write the functional as 




V{w,y) \ dydC,. 



(5.13) 



Let us now perform a one-dimensional monotone decreasing rearrangement of u 
in C for y fixed (see 25 Section II. 3] for details). If u* : S — [0, 1] is the result 
of this rearrangement, then ?2*(-,z) = u{-,z) for all z > R, with R sufficiently 
large, in view of monotonicity of u guaranteed by (5.3) in Proposition |5.4| As 
a consequence, the function u* belongs to the admissible class for . On the 
other hand, by 26 Proposition 12] and [25] Lemma 2.6 and Remark 2.34] we 
know $c[u*] < 'I'clu]: unless u — u* , contradicting the minimizing property of u. 
Then, by strong maximum principle the minimizer u{y, z) is, in fact, a strictly 
decreasing function of z. This proves part (i) of the statement. 

Once the monotonic decrease of u is established, by boundedness of u there 
exists a limit of u{-, z) as z — >■ —oo. The remainder of part (ii) of the statement 
then follows exactly as in ^17^ Proposition 6.6 and Corollary 6.7]. □ 
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Figure 1: A numerical example of non- uniqueness of fronts (see text for details). 



Proof of Theorem 1. The result of parts (i)-(iii) of Theorem [T] follow by combin- 
ing the results of Propositions |3.1[|5.2[|5.4| and |5.5| and taking into account trans- 
lational invariance in z, which allows us to consider only small values of a > 0. 
Also, in (ii) the arguments in the proof of 18, Theorem 3.3] are used to establish 
the estimate on |Vu|. The result of part (iv) follows exactly as in flSl, Theorem 



3.3]. To prove part (v), note that if u satisfies (2.3) and u — w ^ H^{I]), then 
by the assumption on w we have that $™[u] is well-defined, and by minimality 
of $™ [u] we have 

^fiu] = ^:[u] + $f H < + <Pf[w] < ^fiu], (5.14) 

which proves the claim. □ 



6 Global uniqueness of supercritical fronts 



Let us now discuss the question whether the obtained traveling wave solutions 
are unique among all traveling waves with speed c, up to translations. 

First of all, by Proposition 5.4 all bounded positive solutions of (1.31 and 



(1.2) with c > c* must have the asymptotic decay specified by (2.11), which is 
the same decay of the minimizer u of $™. Using the sliding domains method 
of Berestycki and Nirenberg [27], it is possible to show that any two traveling 
front solutions sandwiched between two limiting equilibria must coincide up to 



translation (see 13 Theorem 5.1]). Therefore, in order to have non-uniqueness, 
two different solutions with, say, a = 1 fixed in (3.2), should have different 
limiting behavior at z = — oo. 



Example of non-uniqueness. Clearly, uniqueness of the traveling front 
solutions with prescribed exponential decay cannot hold true in general, even 
in dimension one. Let us demonstrate such non-uniqueness by an explicit ex- 
ample. For simplicity, we present a numerical demonstration of existence of two 
traveling wave solutions with the same value of the velocity c, the result is then 
made rigorous by applying the Theorems of Sec. [2j Namely, let us consider the 
problem (1.3) in the one-dimensional setting with the nonlinear term / given 
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by 



u{l~2u), 0<w<i, 
(2u-l)(l-it)(40u-21), 5<u<l. 



(6.1) 



As can be easily verified, this function satisfies conditions (HI), (H2) and (U). 
One may expect that for a range of velocities c close to cq = 2 the problem ( |1.3[ ) 
with nonlinearity (6.1) has two traveling front solutions: one connecting u — 
and u = ^, and another connecting u = with u = 1. 
for c 



These two solutions 

2.25 are presented in Fig. [ij where Fig. [l](a) shows the traveling wave 
profiles and Fig. [ijb) shows the corresponding phase portraits. From the exact 

traveling wave solution u{z) — 



- J tanh I 



we find that cl 



connecting vi 



with V — I, 



Also, by 
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_ w 2.846 (see Sec. j2jfor the details of the notations). 

Theorem 5.1] we have c* = = 2. Therefore, since cj,^ > c*, 
there are two traveling wave solutions with speed c = 2.25 by Remark |2.2[ one 
connecting zero to vi and the other connecting zero to v. This is what we see 
numerically in Fig. [l] Note that the value of from (2.21) is easily found to 

be c« = w 5.04. 

Let us now present some heuristic arguments as to why the solutions of 
(1.3) and (1.2) are expected to be unique for large enough c. Suppose that 
the nonlinearity / becomes strictly linear for sufficiently small u, i.e. we have 
f{u,y) = fu{0,y)u for all u < 6, with some S > 0. In this case, clearly, the 



decay of any positive solution of (3.1 1 must be precisely 



w{y,z) 



-A_ (c,i'o)z 



■00 (y) 



fc=0 



'0fc(y), 



(6.2) 



for some zq ^ 1- Here we took into account that when < Vk 1^ 0, if ever, 
we have A_(c, z^fc) < A_(c, I'o), hence the presence of a term proportional to 



-A_(c,i/fc)z 



ipk{y) would contradict positivity of the solution. Now, fixing the 



value of a, one can see that any two solutions of the form (6.2 1 differ by a 
function that lies in i7^(E), since all the terms in the series in (6.2) except 
the first one decay exponentially with the rate strictly greater than q- ^ 
the other hand, it is easy to see that for c > c", where is defined in (2.21), 
the functional is strictly convex in iJ^(E). Therefore, fo r ev ery value of a 
and every c > c" there exists a unique positive solution of (1.3 1. This seems 



rather surprising, since, on one hand, uniqueness in this case is a consequence 
of convexity of for c > c", which, in some sense, is dictated by what happens 
in the "core" of the traveling front and thus represents a nonlinear feature of the 
problem. On the other hand, uniqueness also relies on the fact that, regardless 
of the value of c, any two positive solutions with the same asymptotic decay 
have difference that lies in iJ^(S), which is a property of the behavior far ahead 
of the front and represents a linear feature of the problem. From the variational 
point of view, uniqueness of solution is also rather surprising, since, in view of 
Remark|3.4[ there are infinitely many auxiliary functions w that can be used in 
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the construction the functional $Jf' and hence, potentially, different minimizers, 
depending on the choice of w. 

These heuristic arguments turn out to work in the case of generic nonlinear- 
ities in (1.3), with the result stated in Theorem [2j The proof of this theorem 
follows from a sequence of lemmas below. 

Lemma 6.1. We have c" > c*. 

Proof. First, by the definition of c" we trivially have > cq. Now, take any 
u G _ff],(E). Then by hypothesis (H2) there exists w : S -)■ M, with < u < u, 
such that 



|Vu|- - ] dx. 



(6.3) 



Applying the definitions of q and in (2.20), and the Poincare inequality 17 
Lemma 2.1], we obtain 

<i>c«M >l f e^"" {ul + \Vyu\^ - q{y)u') dx 



> 



+ t^o M da; > 0, 



(6.4) 



implying that by (2.101 we have c* < c' 



□ 



Lemma 6.2. Let c> ca and let ui and U2 be two solutions of U.3y and (1.2) 



and whose decay at z = +oo is given by (5.2). Then U2 — Ui € iJ^(S). 
Proof. Introduce h — U2 — ui, which satisfies 

Ah + ch, + f^iO,y)h^if,{0,y)--f,{u,y))h, (6.5) 



for some u bounded between ui and U2. In view of Proposition 5.4 we have 

< u{y,z) < min{l,Ce~^-('=''^°)^}, (6.6) 

for some C > and all {y, z) e E. 

Now, let rj G C^{M.) be a cutoff function with the following properties: ri{z) 
is a non- increasing function of z with super-exponential decay at +00, and 
r]{z) = 1 for z < 0. For fixed e > and c' > let us multiply (6.5) with e'^ '^hrj^, 
where ?7e(z) = rji^z), and integrate over S. After a few integrations by parts, 
we obtain 



\Vh\ 



2 c'(c' - c) 



~ fu{0,y)h^]r,^,dx 



+ j e"^ (2hhzrier]'^ - (c - c)h'^rie'q'^) dx 
e^''{fu{n,y)-fu{0,y))h^Vedx, 



(6.7) 
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where rj'^ = dr]s/dz. 

We are now going to estimate the different terms appearing in (6.7). Let us 
begin by noting that by hypothesis (H2) and (6.6) we have \fu{u,y) — fu{0,y)\ < 
C|u(j/,z)|T < C"(l + eT-^-(=^'^o)^)-i for some C, C" > and all {y,z) G S. 
Therefore, 



On the other hand, since G (E), by Poincare inequality |17, Lemma 2.1] 
it holds 



"-K^l dx > 



e'^'^'-h^rf^ dx 



-2 / e''''hh,'nev'edx~ [ e"' 'h'^\ri'f dx. (6.9) 

Combining this with the previous estimate and using the variational character- 
ization of the lowest eigenvalue vq in (2.6), from (6.7) we obtain 

/2 \ 



cc 



-|- I h^rjl dx 



c\ [ e'''''h^%\r]'^\dx- [ e''''h^\ri'fdx 
Je Je 

<C I + e'^-'^''^''°^')-'^h'^ri^,dx, 



(6.10) 



for some C > independent of e. 

To proceed, we need to choose the cutoff function rj in such a way that the 
terms involving 77^ can be bounded by the terms involving rj^ uniformly in e as 
e — )■ 0. This can be achieved by setting r]{z) = exp(— ^z^) for z > 0. Indeed, 
consider the last term in the left-hand side of ( |6.10| ). We can write 



e''''h^\r)'fdx, (6.11) 



{|i)^|>Qr;} 



for any a > 0. From our choice of r/, the second integral is over all 2; > ae . 
Therefore, assuming h G L^„(I]) for some c" G (0,c'), we have 



e"^/i'|r?;i^c/a; < a' I e^'^h^r^dx 



-2 Jo 



<a^ I e''''-h^rildx + Ce I e''"'h'^dx, (6.12) 
Je je 



where 



Ce = sup 



Ac'-c")zl I 



(6.13) 
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Similarly, 



e'^'^h'^TjeWel dx<a I e^'^h^'^^ dx + a-^Q 



c'zu2 2 , 



e''"'h^ dx. 



(6.14) 



Therefore, we can rewrite (6.10 1 in the following form (also taking (2.91 into 
consideration) 



< Csup 



g(c'-c")z 
1 _|_ g7A_(c, 1^0)2 



+Ce 1 + 



e - c 



e''"^h^ dx 



e''"'h'^dx, 



(6.15) 



for some C > independent of e. 

Now, by an explicit computation, we have 

Ce = sup (a^C^ exp{-e-^{a^C^ - (c' - c")aC))) < a^ (6.16) 
C>i 

as soon as c' — c" < a and e is small enough. Introducing the constants 



(6.17) 



and using (6.16), we then find that (6.151 is equivalent to 



{(c'-c_)(c+-c')-8a|c-c'|} J e'''''h^rildx<C j e''"''h'^dx, (6.18) 

with C > independent of e, whenever c' — c" < a < min(7A_(c, vq), |c' — c|) 
and e is sufficiently small. 

We now note that by Proposition 5.4 we have h G L^„(I]) for some c" G 
(c_,c+). Therefore, we can choose c' — c" + ^ct, with a so small that the 
quantity in the curly brackets in (6.18) is positive. Then, passing to the limit 
e — ^ in (6.18), by monotone convergence theorem we can conclude that, in 
fact, we have h G i^, (E). Applying the same argument to (6.7) we find that 
also h G Hl,{T,), with c' G (c",c+). Finally, iterating this argument, we find 
that 

hGHl,{Y.), Vc'<c+. (6.19) 
In particular, h G i/^(S), which completes the proof. □ 
Lemma 6.3. Let c > cK Then the functional is strictly convex on . 
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Proof. First of all, let us point out that due to the presence of the exponential 
weight the functional ^™ is not twice continuously differentiable in H^(T,), so 
the convexity of the functional cannot be concluded by evaluating the second 
variation of Nevertheless, what we will show below is that the functional 
Vl/™ is strongly convex 29 , i.e., there exists a > such that for every /ii, /12 € ■A.'^ 
it holds 

^:[h2] > ^:[hi] + 5^:{h,M - hi) + |||/i2 - /iilli.(s)- (6.20) 
Indeed, by hypothesis (H2) we have 

+ W e''[^y{h2-hr)\^~f^(Si,y){h2~h{f)dx, (6.21) 



for some < u < 1. Now, by Poincare inequality [17, Lemma 2.1], we have 

^2 



e^ia,(/i2-/;,i)|- 



dx> — Who 
- 8 " 



and by the definition oivQ in (2.201 



(|V,(/i2 



i|2 



^0 I 



fuiu,y){h2 - hi)"^) dx > y 11/^2 



L?(S) 



(6.22) 



(6.23) 



1 ^2 



> 0, in view of c > c" and 

□ 



Strong convexity then follows with a = 
implies strict convexity 29 Chapter 12. H*]. 

Proof of Theorem 3. Let ui and U2 be as in Lemm a|6.2[ Without loss of gener- 
ality we may assume that a is sufficiently small in (5.2). Let us define w to be a 
suitably truncated version of u i. Then, by hypothesis (H2) the functional ^E"^ is 
of class and by Lemma 6.3 is also strictly convex on . Then hi := ui — w 
belongs to . In fact, the bmction /12 ■= U2 — w also belongs to , in view 
of Lemma 6.2 On the other hand, since the equivalent functional <I>^ is also 
strictly convex, it has at most one critical point. Therefore, since ui and U2 solve 
the Euler-Lagrange equation (1.3) for $™ and are, therefore, critical points of 
we conclude that ui = M2. □ 



7 Global asymptotic stability of fast fronts 

We now investigate the stability properties of traveling wave solutions with 
speed c > c* obtained in the preceding sections. Let us point out in the first 
place that one cannot expect these traveling wave solutions to be stable in 
the usual sense (say, with respect to small L°° perturbations). Indeed, under 
hypotheses (HI), (112) and (U) a truncation of the traveling wave solution with 
any c > c* by a cutoff function that vanishes identically for all z > i? 3> 1 would 



27 



not converge to the corresponding traveling wave, but will instead propagate 



with the asymptotic speed c* of the critical front 18 Theorem 5.8]. Similarly, 
modifying the rate of exponential decay in front of the traveling wave solution 
with speed c > c* may lead to acceleration or slowdown of the front (see e.g. 
30 -32]), or even lead to irregular behavior 33 . Therefore, stability of these 



traveling wave solutions should be considered with respect to suitable classes of 
perturbations 30p4 . Here we show global stability of the (unique) supercritical 
traveling waves with c > c" under perturbations lying in the weighted space 
L^, (E) with suitably chosen c' > 0. 



We begin with a general existence result for the solutions of ( 1.1 ) of the form 



u{y, z, t) ~ u{y, z — ct) + w{y, z — ct, t), 



(7.1) 



where u is given by Theorem [T] with some fixed c > c* , and w lying in an 
appropriate exponentially weighted space. The function w satisfies the following 
parabolic equation 



wt = Aw + cwz + !{u + w,y) - f{u, y). 



(7.2) 



Easily adapting 20 Proposition 3.1] (see also |18|), we can obtain the following 
basic existence and regularity result for the initial value problem associated with 



(7.2): 



Proposition 7.1. For any c' > 0, let wq G ic'(^) '^'^'^ let < u + w < 1. Then 
there exists a unique solution w £ C"((0, oo); iJ^, (S)) n C^'"((0, +oo); L2,(E)) n 
C{[0,+oo);L'^,{'S)), for alia £ (0,1), of ( |7.2[ ) and the same boundary conditions 
as in ( |1.2[ ), such that w{x,0) — wo{x). 

We are now in a position to complete the proof of Theorem [3j following the 
method introduced in [Tg] . 

Proof of Theorem 3. Without loss of generality, we may assume that c' < 2A+(c, i>o), 
where A±(c, z^) are defined in (2.11). The result in Proposition 7.1 implies that 
we can multiply (|7.2| hy e'^ '^w and integrate over S. Using hypothesis (H2) and 



the definition of q{y) in (2.20), we have 



e'^^wwtdx< / e'^ ^ (Aw + cw z — q{y)w) w dx . 
Then, after a number of integrations by parts, we obtain 



(7.3) 



'w^ dx<-2 e^'lwl 



c'(c-c') 



+ \yyw\' ~ q{y)w' dx. (7.4) 



Applying the Poincare inequality fl7' Lemma 2.1] and the definition of i)o in 



(|2.20|), we then arrive at 

c'z„..2 



dtj^ - I 4 2 



1 



{c'-2A_(c,;>o)}{c'-2A+(c,i>o)} / e'^w^da; 




e'^'^w^ dx 



(7.5) 
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Therefore, since c' S (2A_(c, i>o), 2A+(c, ;Jo))i the coefRcient in front of the inte- 
gral in the right-hand side of (7.5) is negative, which impUes exponential decay 
of the iy^,-norm of w. □ 
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